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Abstract. We consider the possibility that the discrete long-range ordered states of Er2Ti207 
are selected energetically at the mean field level as an alternative scenario that suggests selection 
via thermal fluctuations. We show that nearest neighbour exchange interactions alone are not 

sufficient for this purpose, but that anisotropics arising from excited single ion crystal field 
states in Er2Ti207, together with appropriate anisotropic exchange interactions, can produce 
the required long range order. However, the effect of the single ion anisotropics is rather weak 
so we expect thermal or quantum fluctuations, in some guise, to be ultimately important in this 
material. We reproduce recent experimental results for the variation of magnetic Bragg peak 
intensities as a function of magnetic field. 



1. Introduction 

Frustration is a common feature of magnetic interactions on a pyrochlore lattice. It often has the 
consequence of introducing a large degree of degeneracy into the ground state. This might be 
relieved by weaker interactions with the result that long-range order sets in at temperatures much 
smaller than the overall energy scale set by the interactions which, itself, is typically correlated 
with the Curie- Weiss temperature. Fluctuations alone can relieve the energetic degeneracy (e.g. 
[1, 2]) by a mechanism called order-by-disorder or even select states that do not lie within the 
ground state manifold [3, 4]. The suppression of the transition temperature compared to the 
characteristic energy scale of the interactions can result in quantum fluctuations being more 
important than thermal fluctuations in dictating the ordering. Another possibility is that the 
degeneracy is not relieved in spite of weaker interactions or fluctuations so that spin hquids or 
states with unconventional long-range order emerge (see, for example [5, 6, 7]). 

The diversity of theoretical possibilities that arises as a consequence of magnetic frustration is 
paralleled by experiments which have brought many unusual materials to light. Even among the 
heavy rare earth titanate pyrochlore magnets, one can find materials that have no conventional 
long-range order at the lowest temperatures explored - the spin ices Ho2Ti207 and Dy2Ti207 
[8], and the cooperative paramagnets Tb2Ti207 [9] and Yb2Ti207 [10], and materials which do 
exhibit conventional long-range order, but which are puzzling in their own right: Gd2Ti207 [11] 
and Er2Ti207 [12, 13, 14]. 



Figure 1. One ordered domain {ip2 state) 
of Er2Ti207. The spins are in the easy 
planes perpendicular to the local (111) 
directions. The other five possible ground 
state configurations are obtained by discrete 
rotations and time reversal. 

In this article, we consider the latter material. Er2Ti207, which has a Curie-Weiss 
temperature ^cw ~ —15.9 K [15], undergoes a zero-field transition at about 1.2 K into a 
phase with ordering wavevector q = [12, 13]. The magnetic structure has been characterized 
through an analysis of spherical neutron polarimetry data [13]. There are six domains which can 
be represented by showing the spin configuration on a single tetrahedron. Figure 1 represents 
the spin orientation on a tetrahedral primitive unit cell for a q = ordered state; two others 
are obtained by discrete global rotations and the remainder are time reverses of the first three. 
Henceforth, we refer to these as 1^2 states to conform to the notation used in [12] and [13]; 
these states transform according to the magnetic two dimensional irreducible representation 
Eg based on the octahedral point group Oh of the pyrochlore lattice for which the spins are 
oriented perpendicular to the local (111) directions. The crystal field is responsible for the local 
xy anisotropy, but the occurrence of discrete ordered states is not presently understood [2, 12]. 

An inelastic neutron scattering scan for scattering along the q = [22L] direction has revealed 
a gapless or almost gapless mode in zero field at q = [220] [14]. As a magnetic field H along 
[110] increases from zero, one finds that the spectrum is gapped for fields around |H| ~ 0.5T and 
for strong fields |H| ~ 5.0T but that there is a regime in between where it is apparently gapless. 
It has been suggested, on this basis that there is a field-driven quantum phase transition in this 
material [14]. 

In connection with Er2Ti207, an easy plane model, with classical spins constrained to lie in the 
local (111) planes and interacting via isotropic exchange [2, 12], has been studied. Energetically, 
this model has an extensive ground state degeneracy but Monte Carlo simulations show that 
there is a thermally-driven order-by-disorder transition to the ^^2 states. However, it is known 
that this cannot be the complete solution to the problem of 1^2 selection in Er2Ti207 [2, 12]. 
First of all, experimentally, the zero field transition is known to be continuous, whereas the 
simulations reveal a strong first order transition. Another problem is that the Er^"*" ions have 
a dipole moment that breaks the degeneracy associated with the isotropic exchange, leading to 
a different q = ordering to that seen in Er2Ti207 [16]. Finally, the heat capacity in the low 
temperature phase varies as whereas the model predicts a linear variation with temperature. 

This paper investigates the conditions under which the selection of 1^2 states as ground states 
may be determined energetically and at the mean field level in the real compound. We reason 
as follows: regardless of the mechanism behind the observed ordering, a six- fold anisotropy 
must be present in the real system to select the discrete set of states. Whether this anisotropy 
is generated "energetically", "entropically" or "dynamically" via quantum fluctuations is most 
likely not a meaningful question from an experimental point of view because, in all likelihood, the 
experimental observables determined for a material exhibiting true order-by-disorder could be 
mimicked by an appropriate choice of interactions that uniquely select the ordered states. Given 
that the only known order-by-disorder mechanism is problematic and that energetic selection can 
be studied, for example, within mean field theory, we adopt the latter approach. The question 
is then whether such interactions could be natural in Er2Ti207 in the sense that they could be 
justified within a microscopic model. We investigate this issue in the following. 




2. Selection by exchange interactions 

In this section, we look at the problem of whether the most general nearest neighbour 
interactions, expressed in terms of the Er'^+ angular momenta on a pyrochlore, can break the 
local XY spin rotation symmetry down to the observed ordered states of Er2Ti207 without any 
single ion anisotropies. We show that, in fact, these interactions are not sufficient to select 
the required states. We start from the nearest neighbour exchange interactions allowed by 
symmetry on a pyrochlore lattice which can be found in [17] and which we present here, in 
Table 1, with respect to four sets of local axes - one for each sublattice. The sum of the terms 
in each of the four columns of the Table is a linear combination of bilinear operators that is 
invariant under the point group of the pyrochlore lattice Oh- We denote the four invariants 
by Xp for p = 1,...,4. The isotropic exchange interaction, nearest neighbour dipole-dipole 
interaction, Jj ■ Jj — 3(Ji • rij){Jj ■ Vij), where f are unit vectors connecting nearest neighbours, 
and the Dzyaloshinskii-Moriya (DM) interaction, D • ( Jj x Jj) [18, 19], can be written as linear 

combinations of these invariants. For example, isotropic exchange J2{ij)'^i ' = J2p=i-^p- 
Ordinarily the dipole-dipole interaction is long-ranged; because we cut it off beyond nearest 
neighbours, which is not unreasonable for modelling a zero moment q = ordering, we refer to 
it as the pseudo-dipole interaction. 

We can express these exchange interactions in terms of bilinears that are the products of 
basis vectors of the point group. For Er2Ti207, the irreducible representation that goes critical 
at the transition temperature is the two dimensional irreducible representation Eg [12] which 
has orthogonal basis vectors that we denote by Je^ and Je_ from which the order parameter 
can be constructed. These basis vectors are written explicitly in terms of J" components in 
equation A.l in the Appendix. One can show that these basis vectors appear in the symmetry 
allowed exchange interactions Xp only in the combination Je^ Je^ (see Appendix) . The crucial 
point of the argument is that this exchange term does not distinguish any directions in the easy 
planes of the four sublattices so, in the thermodynamic limit, one would expect a continuous 
symmetry breaking provided the coefficient of Je+Je^ in the Hamiltonian is negative. This 
situation is captured by a Landau theory based on the Eg irreducible representation to quartic 
order in the order parameter. The order parameter for Er2Ti207 is the expectation value of 
ot E+ Je+ Je- for = (XE- = ^ and then five other values related to the first by symmetry. 
These are minima of a Landau theory to sixth order in the invariants of Eg (see, for example, 
[20]). 

The conclusion of this section is that nearest neighbour exchange interactions allowed by 
symmetry on a pyrochlore cannot, in the absence of single ion anisotropies, break the symmetry 
uniquely down to the -02 states. In the next section, we explore a mechanism for the generation 
of sixth order terms in the Landau theory for Er2Ti207. 

3. Inclusion of in-plane anisotropies 

We mentioned in Section 1 that, with isotropic exchange interactions, the pyrochlore with 
XY anisotropy has an extensive ground state degeneracy [2], and that the discrete ip2 states 
belong to this continuous set of ground states. Then, the inclusion of a six-fold anisotropy, 
with minima along the spin directions of the ^/'2 states, must break the degeneracy exclusively 
down to the six tp2 domains observed in Er2Ti207. Single ion anisotropies were omitted in 
the symmetry arguments of Section 2. With classical spins, the relevant single ion anisotropy 
is + ~ cos(6^), where ^ is an angle in the local (111) plane. This is proportional to 
the Racah operator Og which appears in the crystal field Hamiltonian [21]. So, although the 
lowest crystal field doublet has almost perfect XY degeneracy, one might expect that the mixing 
of excited crystal field levels into the ground state induced by the interactions between the 
angular momenta could perhaps introduce a six-fold anisotropy of the correct modulation to 
break the symmetry energetically. This would fulfil the requirement that the Hamiltonian be 



Table 1. Exchange interactions over a single tetrahedron expressed in terms of local coordinates 
for each magnetic ion. The coordinate system is described in the Appendix. Invariant Xp (for 
each p) is the sum over terms in its column. Term Jj ■ 3j for given i and j is the sum over terms 
in its row. e = exp(27ri/3). h.c. stands for Hermitian conjugate. 
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well-motivated from a microscopic point of view. We note, however, that recent inelastic neutron 
scattering experiments on Er2Ti207 [14] do not resolve a gap to spin excitations which should 
be present if the symmetry is broken in the ground state. However, as we point out below, one 
would expect this gap, if it exists, to be small. 

To investigate whether the excited crystal field levels can introduce a six-fold anisotropy with 
the correct modulation, we have studied a microscopic Hamiltonian H = Hd + Hex where 

^^cf = E ^oOo. + B'oOl, + i?|0|, + B^Ol + Bpl + bIoI (1) 

i 

^ex = ^ E E E Jafi^,M^,)J^{^.)- (2) 

/i,!^ a,b a,j3 

The Greek indices are angular momentum components, a and h are sublattice indices (ranging 
from 1 to 4) and is a primitive lattice vector with Yi^^ = — Rj^; index i runs over all 
magnetic ions. The crystal field parameters are rescaled from those of Ho2Ti207 [21], and the 
choice of nearest neighbour exchange couplings iOp is guided by the requirement that, in the phase 
diagrams to be presented below, the order of magnitude of the couplings should be consistent 
with the experimental Curie-Weiss temperature [15] which we have fitted by computing the 
uniform susceptibility in the random phase approximation (RPA) [22], and in the paramagnetic 
phase, from the Hamiltonian given above. We rewrite ilex on a single tetrahedron as Yl'^=i '^p^p^ 

hence constraining the elements of J^^ ■ 

To proceed, we carry out a mean-field decoupling of J^{ix)J^{u) by writing it as 

and omitting the last term which quantifies the fluctuations. The angled brackets indicate a 
thermal average (^4) = -|.Tr \^Ae~^^^, and we have abbreviated J^(R^) by J^ilA- 
The original interacting problem is reduced to a single ion problem 

IJ,,a,a ^ ' 

where is an effective magnetic field of the form ^"(/i) = Ylub i3'^af(^i^'')('^b('^))- 



Mean field solutions for a cubic unit cell and examination of the RPA susceptibility show that, 
for the parameters we have explored, the ordering wavevector is q = 0. So, for a qualitative 
investigation, we can, at this time, restrict our attention to a single tetrahedron. 

Starting from random choices for the twelve expectation values (J^), the Hamiltonian is 
iterated to self-consistency for different couplings. We have found that the six ?/)2 domains 
arc solutions to the mean field equations for some regions of parameter space. Figures 2 and 
3 show phase diagrams including each showing such a region - isotropic exchange and pseudo- 
dipole interactions, and isotropic exchange with DM interactions respectively. The self-consistent 
solutions are indicated on the phase diagrams. As expected, we find the q = states reported in 
[16] with isotropic exchange and positive pseudo-dipole coupling. Though the physical dipole- 
dipole coupling is positive, we find the Er2Ti207 ordered states only when the pseudo-dipole 
coupling is negative in this case. We find the same ^"2 states with positive DM coupling and 
isotropic exchange for J < Jj)M (see Fig. 3). We have also found couplings xi > and X2 < 0, 
given in the Appendix, such that V'2 order occurs for the positive pseudo-dipole coupling of 
D = 0.02 K appropriate to Er2Ti207[12]. 



The progress of iterations of the mean field equations to self-consistency sheds some light 
on the robustness of the energetic selection mechanism. Figure 4 shows the free energy per 
site plotted against the number of iterations when the couplings are such that the ip2 states are 
selected. The first few iterations have been cut off to see the scale of the trend. An accompanying 
plot. Figure 5, shows the variation in the angle of {J) in the local (111) plane. The spins settle 
into their respective xy planes within a few iterations but around O(IO^) iterations are required 
for the spins to converge to the ^^2 angles. The variation in the free energy across this interval 
is 0{W~^) K. This is in marked contrast to the case of isotropic exchange and positive pseudo- 
dipole coupling where the convergence is complete within O(IO^) iterations corresponding to 
energetic selection without any need for single ion anisotropics within the local xy planes coming 
from excited crystal field levels. 

This result can be understood in terms of a mean field treatment of the quantum corrections to 
the Hamiltonian obtained by projecting onto the space spanned by the ground state crystal field 
doublet (see, for instance, [23]). Such corrections are computed perturbatively in powers of J/ A 




Figure 2. Phase diagram for isotropic 
exchange J and pseudo-dipole D couplings. 
Three different q = phases arise as 
indicated. From bottom to top, there a 
configuration with a net moment, then the ■02 
state, and then the state reported in [16]. 



Figure 3. Phase diagram for isotropic ex- 
change J and Dzyaloshinskii-Moriya interac- 
tions Jdm- From bottom to top, there is a 
state familiar for dipole interactions [16], then 
the '02 state and then a ip3 state [12]. 



where J is the characteristic energy scale of the exchange, which is about 0.1 K for Er2Ti207, 
and A ~ 10^ K is the scale of the crystal field splitting [12]. These terms involve admixing 
with excited crystal field states which are responsible for bringing the six-fold anisotropy into 
the mean field solutions. The smallness of J/ A and the fact that six- fold anisotropy cannot be 
generated to lowest order in the perturbative corrections implies that the selection of ijj2 states 
is very weak within the mechanism considered here. 
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Figure 4. Variation of the free energy per 
tetrahedron upon iteration with J = 0.1 K 
and D = -0.02 K. T = 0.3 K. 
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Figure 5. Convergence of the angle upon 
iteration with mean field parameters J = 0.1 
K and L> = -0.02 K. T = 0.3 K. 



4. Bragg scattering in a magnetic field 

We have studied the effect on the elastic neutron scattering pattern of a magnetic field applied to 
the material in its low temperature ordered state. Following experiments by Ruff and co-workers 
[14] in which the intensities of five Bragg peaks were measured as the field strength was varied, we 
have applied a field in the crystallographic [110] direction. A Zeeman term Hz = —fiBd'^'H was 
added to .ffMF (Equation 4). The field was increased incrementally starting from one zero field 
1^2 domain and the self-consistent solution was found for each field. Figure 6 is the mean field 
result and Figure 7 are the experimental intensities taken from [14] . In the experimental range 
|H| = T to |H| = 5 T, there is good qualitative agreement between the two sets of intensities 
for the chosen couplings. The maximum |H| = 5 T is sufficiently small that the spins remain 
in their respective xy planes. The change in the spin configuration coincides with the picture 
presented in [14] that has spins along chains in the field direction maintaining antiferromagnetic 
correlations before rotating together into the field direction. Most of the spin rotation happens 
in a narrow band around 1.5 T in both the experiment and in the theory. However, one can 
change the band of fields for which this sharp feature is seen by varying the couplings (of course, 
with the condition that the ^^2 states are the H = ordered states). 

5. Discussion 

We have asked in this article whether models for Er2Ti207 can be written down that are 
microscopically well-motivated and which energetically select the experimentally observed ip2 
states. This approach stands in contrast to earlier investigations that focus on fluctuation- 
induced ordering [2, 12]. We have found that there are indeed such models with crystal 
field terms and suitable symmetry-allowed bilinear interactions between the nearest neighbour 
angular momenta. We have shown that models with only bilinear interactions cannot select the 
required states. The same is true for local (111) plane xy models with these bilinear interactions. 
The introduction of crystal field terms does allow for 1^2 states to be selected for certain regions 
in the space of couplings - the interaction-induced admixing of the ground doublet with excited 
crystal field levels being the source of the necessary six-fold anisotropy. However, because the 
energy scale of the crystal field splitting is very much greater than the characteristic scale of the 
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Figure 6. Bragg peak intensities from mean 
field theory with isotropic exchange J = 0.1 
K and dipole-dipole couphng D = —0.02 K. 
T = 0.3 K. 



Figure 7. Experimental magnetic Bragg 
peak intensities for Er2Ti207 in a [1,1,0] 
field taken from [14]. 



interactions, the selection of these discrete ground states is extremely weak. Even if the exchange 
interactions were such that 1^2 long-range ordering would occur, one would expect the transition 
temperature Tc to be very much smaller than the experimental Tc provided that an order-by- 
disorder mechanism does not select the ■02 states. In view of these results, and the problems 
with the isotropic exchange order-by-disorder mechanism that were outlined in the Introduction, 
we tentatively conclude that neither of these proposals is sufficient on its own. It is conceivable, 
instead, that they complement one another - that the microscopic model might weakly select the 
1^2 states and that fluctuations have the effect of raising the transition temperature. One may 
ask whether order-by-disorder could operate with a continuous transition rather than a strongly 
first order transition in these circumstances. It might also be that quantum fluctuations must 
be considered to find the correct answer to this question as suggested in [12]. 

One success of our mean field theory approach is that we have been able to reproduce the 
experimental variation of the magnetic Bragg peaks intensities in a magnetic field in the [110] 
direction. Provided we choose the interactions correctly, in the range of fields where the Bragg 
intensities change most sharply, the theory matches the experimental field strength. The theory 
does not reproduce the upturn in intensities between T and 0.5 T however; the mean field 
theory selects two domains out of the six in very weak fields (|H| <C 0.5 T). Two interpretations 
for this upturn have been suggested in the literature: the first is that it is due to the selection 
of preferred domains [12]. The second interpretation comes from the observation that a broad 
component to the scattering is suppressed for weak fields and that this coincides with the upturn 
in the main scattering component to conserve the intensity. It has been suggested that quantum 
fluctuations are responsible for the broader scattering peaks and that these are suppressed 
rapidly in weak fields [14]. Such an effect may lend weight to an investigation into zero field 
ordering in Er2Ti207 induced by quantum effects. 
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Appendix 

The basis of four magnetic ions with respect to face centred cubic system is ri = |(0,0,0), 
r2 = 1(1,1,0), r3 = 1(1,0,1) and r4 = |(0, 1,1) where a is the edge length of the cubic unit 
cell. The local coordinate system is chosen with local z axes in the (111) directions. The x 
axes are chosen as follows xi = ;^(1,1,— 2), X2 = ;^(— 1, — 1, — 2), X3 = ;^(— 1,1,2) and 

X4 = ;^(1, — 1, 2) and the axes are right-handed. 

The basis vectors for the different irreducible representations of are given below. The four 



local z components of the angular momenta are decomposed into T\ and the remaining 



components are E ®T\®T2. 

Je+ = Ji + J2 +J3+J4 Je_ = Ji + J2 + -^3" + J4 (A.l) 

Ja, = JI + J| + J| + J| (A.2) 

■^Ti.i = Ji ~ Ji ~ J4 Jti^i = Ji — Ji J3 ~ (^-3) 

j^,^, = J! + Ji-Ji-J! ' (A.4) 

= - 4 - 4 + 4) + h-c. = \e{Jt - Jt + 4 - 4) + h.c. (A.5) 

Jk,, = 4 + J2-4-4 (A.6) 

= -l^*i4 -4- 4 + 4) + h-c. 4, = - Jt + Jt - 4) + h-c. (A.7) 

Jk = 4 + 4 - 4 - 4 (A.8) 



The appearances of the Ti irreducible representation have labels to distinguish them: Ti^i and 
Ti^2- The phase e = exp(27ri/3). The exchange invariants, given in Table 1, can be written in 
terms of these basis vectors as follows 

1 2 1 2 V2 
^1 = + ^•'n,! ^2 = — §"''^1.1 ■ ''^1,2 

^3 = \^T,,2 - l^T2 ^4 = -^Je+JE- + ^Jri,2 + ^Jt2 

SO, as stated in Section 2, the only appearance of the basis vectors of the two dimensional 
irreducible representation Eg, which control the symmetry breaking to the 1^2 states, is in the 
combination Je+Je- which does not lead uniquely to the V'2 states, as outlined in Section 2. 

In Section 3, we make reference to invariants {x\, which are linear combinations of those in 
Table 1, xi = 2Xi + \X2 + 5X3 + 2X4 and X2 = -Xi + \X2 + ^X^ - X4, which have the 
property that xi + X2 is the isotropic exchange on a tetrahedron. 
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